Kosterlitz-Thouless phase transition and reentrance in an anisotropic 3-state Potts 

model on the generalized Kagome lattice 
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The unusual reentrant phenomenon is observed in the anisotropic 3-state Potts model on a gen- 
eralized Kagome lattice. By employing the linearized tensor renormalization group method, we find 
that the reentrance can appear in the region not only under a partial ordered phase as commonly 
known but also a phase without a local order parameter, which is uncovered to fall into the uni- 
versality of the Kosterlitz-Thouless (KT) type. The region of the reentrance depends strongly on 
the ratios of the next nearest couplings a — J2/IJ1I and /3 = Js/|Ji|. The phase diagrams in the 
plane of temperature versus /3 for different a are obtained. Through massive calculations, it is also 
revealed that the quasi-entanglement entropy can be used to accurately detect the KT transition 
temperature. 

PACS numbers: 75.10.Jm, 75.40.Cx 



I. INTRODUCTION 

Two-dimensional (2D) frustrated classical lattices pos- 
sess many interesting physical properties (see, e.g. 
Refs. [1—8]. One of them is the so-called reentrant 
phenomenon 1,3,6,9 , which is defined as the occurrence of 
a disordered phase (usually a paramagnetic phase) in the 
region under an ordered or a partial ordered phase on the 
temperature scale. This disordered phase is called the 
reentrant phase that has been detected experimentally in 
spin glasses 10 and studied theoretically in several exactly 
solvable 2D frustrated Ising systems on lattices such as 
the centered square lattice 1 , generalized Kagome lattice 3 , 
centered honeycomb lattice 2 and other three-dimensional 
lattices 5 . Up to now, the reason for the occurrence of the 
reentrant phase is still under debated. 

Some conjectures 1 ' 3,5 have been proposed to under- 
stand the reentrant phenomenon, in which the essential 
is that the reentrant phase is probably caused by frus- 
trations and in the ground state, there should be at least 
one partial order phase adjacent to an ordered phase or 
another partial ordered phase between which the reen- 
trant phase can appear at finite temperature. Above 
the reentrant phase, there is usually a partial ordered 
phase whose disordered sublattices supplement the en- 
tropy which is lost due to the formation of the ordered 
sublattice. In addition, other factors such as the coordi- 
nate numbers 5 of the sites on the disordered sublattice 
and the freedom of the on-site spin 9 (e.g. the value of q in 
the Potts model) may also have effects on the reentrant 
phenomenon. 

Alternatively, to explore the incentives of the reen- 
trant phenomenon, the q-state Potts model 11 can offer 
some clues. Potts model is the generalization of Ising 
model by extending the on-site freedom from q = 2 to 
q > 2. The q-state Potts model does not have exact so- 
lutions. For some lattices such as the piled-up domino 
model (frustrated Villain lattice) 9 , the Potts model pos- 
sesses a reentrant phase that could maintain in the region 



under one partial ordered phase when 1 < q < 4 (except 
for q = 2), which indicates that the occurrence of reen- 
trant phenomenon has close connection to the value of 

q- 

In this work, we shall focus on the anisotropic 3-statc 
Potts model on a generalized Kagome lattice, as shown 
in Fig. 1, which contains three different couplings: J\ 
is the diagonal coupling and is presumed to be ferro- 
magnetic (F), J2 and J3 are couplings along the verti- 
cal and horizontal directions, respectively, both of which 
can be either antiferromagnetic (AF) or F. For the Ising 
case {q = 2), in the ground states there exist two par- 
tial disordered configurations A and B, as shown in Fig. 
1. Configuration A corresponds to the case of J3 < 
(F) and J 2 > (AF), and B is the case of J 2 , J3>0. 
The spins on the central site in both cases are in the 
free states, therefore the order is defined as partial order 
or partial disorder. By the exact solution the reentrant 
phase is found between the F phase and the partial or- 
dered phase with configuration A 3 , which is induced by 
frustrations. It is interesting to ask whether there is any 
reentrant phenomenon in the present anisotropic 3-state 
Potts model on the generalized Kagome lattice, in which 
J2 is set to be negative and J3 positive with J\ = -1 (F). 
Such a mixed (with AF and F interactions) Potts model 
cannot be solved exactly, and we perform numerical sim- 
ulations to obtain the specific heat, susceptibility, cor- 
relation length, and entanglement entropy of the system 
to determine the phase diagram precisely It is discov- 
ered that although the thermal fluctuation depresses any 
partial order with configuration A and B in the ground 
state, our numerical results strongly support that there 
still exists a reentrance but with the Kosterlitz-Thouless 
(KT) type phase transition at certain values of Ja/|Ji| 
and Jz/\Ji\- In addition, it is found that the reentrance 
disappears when q > 3. 

This paper is organized as follows. In Sec. II, the 
model Hamiltonian and the partition function are de- 
fined, and the tensor network representation of this 
model is given. In Sec. Ill, the specific heat, susccpti- 
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bility and the entanglement entropy are calculated, and 
the KT type phase transition is discussed. In Sec. IV, by 
examining the singularities of the thermodynamic quan- 
tities, the phase diagrams in (3 — T plane are presented 
for a = — 1, < a < 1 and a > 1, respectively, and the 
reentrant phenomena are observed. Finally, a summary 
is given. 

II. MODEL AND PARTITION FUNCTION 

The Hamiltonian of the anisotropic q-state Potts 
model on the generalized Kagome lattice reads 

(»>j) 

where Ji (i = 1, 2, 3) are coupling constants, S ai(Tj is the 
Kronecker symbol with cr.; = 1, 2, • • • q. In the following 
we shall focus on q=3. For clarity, we define a = Jij\J\ |, 
and/3 = J 3 /|Ji|. 

By means of the Trotter-Suzuki decomposition the par- 
tition function of the system can be written as 

Z = Tre- H ' T = Tr\[e-* h : (2) 
h = Ji(8 aai + S aa2 + 5 aa3 + 5 a<Ji ) (3) 

+M5 ) + + <WJi 

where T is temperature, the Boltzmann's constant fc^ is 
taken as unity, and e is the Trotter slice. 

To study this present anisotropic 3-state Potts model, 
we employ the linearized tensor renormalization group 
(LTRG) method 12 , a recently developed numerical al- 
gorithm for calculating the thermodynamic properties 
of the low-dimensional quantum lattice systems with 
highly accuracy and efficiency, which has been success- 
fully applied to a few one-dimensional quantum spin lat- 
tice models 13,14 . In the tensor network representation 
of the partition function (Fig. 2), e~ th can be seen as a 
fourth-order tensor T as shown in Fig. 2(b), and the par- 
tition function Z can be represented as a tensor network 
in Fig. 2(c). During the LTRG calculations, we keep the 
bond dimension cut-off D c of the tensor network at least 
60 and the truncation error is less than 10 -7 . 



III. THERMODYNAMIC QUANTITIES AND 
KOSTERLITZ-THOULESS PHASE TRANSITION 

Kosterlitz-Thoulcss phase transition in the mixed 
Potts model 15,16 is characterized by the divergence of the 
correlation length £ near the critical temperature T Ci kt 
of the form 

£ ~ e v""c!w j (4) 

while the specific heat shows no divergence and has a 
broad peak near T Ci kt- Below T^^t, there is a phase 




FIG. 1: (Color online) The geometric structure (left) and 
partial order configurations (right) of the generalized Kagome 
lattice. The thin blue line represents ferromagnetic coupling 
Ji, the thick solid line (J3) and slash black line (J2) could be 
either ferromagnetic or antiferromagnetic coupling. The black 
and red dots in partial order A and B configurations indicate 
different directions of the on-site spin and the centered green 
site is a free spin. 

(coined as floating phase) with an algebraically decaying 
correlation function. Such kind of phase was explained 
as the resemble of vortices that are closely bound in pairs 
in 2D XY models 16 . In the present 3-state Potts model, 
we also found such a phase and in the following we will 
follow Ref. [15] to call it the floating phase. 

Now let us first explore the thermodynamic properties 
of the generalized Kagome Potts model with parameters 
a = J2/IJ1I = -1 and (3 = Jz/\J\\ > 0. The tempera- 
ture dependence of the specific heat C for different /? is 
obtained, as shown in Fig. 3. It is seen that, for j3 < 
1, the specific heat has sharp peaks at critical temper- 
atures [Fig. 3(a)], which indicate that the second-order 
phase transitions may take place in this case. This could 
be easily understood that when (3 < 1 the ferromagnetic 
couplings J2 and J\ are dominant, and for a given /3 < 1 
the system may undergo a phase transition from ferro- 
magnetic phase to paramagnetic phase, which will also be 
confirmed in calculations of magnetization (see Fig. 4). 
It is observed that the critical temperature T c decreases 
with increasing f3 when /3 < 1. For > 1, the spe- 
cific heat shows broad peaks at low temperatures, and no 
divergence is found, displaying that no thermodynamic 
phase transitions can happen in this case. However, this 
does not imply that the topological phase transition like 
KT phase transition is unlikely. In fact, it is the case, as 
seen below. 

To understand further the thermodynamic behavior of 
this model, we come to look at the magnetization per site 
defined as 

i 

where N is the number of lattice sites, and (...) denotes 
the thermal average. Note that in this definition m de- 



scribes the mean value of magnetization on each site 
whose value equals | when the lattice is in disordered 
state. Figure 4 presents the temperature dependence of 
magnetization per site at a = — 1 for different (3 < 1 . The 
sharp changes of m at critical temperatures can be seen, 
showing that for a given (3 < 1 the system indeed has a 
thermodynamic order-disorder phase transition with in- 
creasing temperature, being consistent with the results 
of the specific heat [Fig. 3(a)]. For f3 > 1, m = at fi- 
nite temperature, showing that in this case no long-range 
order appears. To examine if there is a topological phase 
transition in the case of (3 > 1, we calculated the suscep- 
tibility x °f the system defined by 



dm 
~dh' 



(0) 



where h is the uniform magnetic field. The results are 
given in Fig. 5, where, for different values of j3 > 1, there 
are singularities (kinks) at certain temperatures, imply- 
ing that there might be a kind of non-thcrmodynamic 
phase transition occurs. 

To confirm such topological phase transition is of the 
KT type, we investigate whether the correlation length 
£ has the form of Eq. (4). Following the lines intro- 
duced in Ref. 17, we obtained the results of correlation 
length versus temperature, as shown in Fig. 6 (a) . Ap- 
parently, when (3 > 1 near the transition temperature 
T C .KT, the correlation length £ shows the same behavior 
as that given in Eq. (4), where a slight deviation could 
be modified by increasing the bond dimension D c . This 
result supports that the phase transition for /3 > 1 is in- 
deed of the KT type. For j3 < 1, the correlation length 
£ reveals a distinct behavior from the case of (3 > 1, as 
shown in the inset of Fig. 6 (a), consistent with the pre- 
vious observation that the phase transition occurring in 
the case of j3 < 1 is thermodynamic. Besides, by fitting 
the curve, we obtain the critical exponent v = 0.820 for 
(3 < 1 in this model, which is close to the conjectured 
value | (Ref. [11]). 

In addition, by making use of the elements in the A ma- 
trix in the matrix product state (MPS), shown in Fig. 2 
(d) , obtained by the contraction of the tensor network in 
Fig. 2 (c), we can introduce the quasi-entanglement en- 
tropy following the similar way in one-dimensional quan- 
tum lattice system 



Sq = -T, i X i logX l . 



(7) 



Actually, the MPS obtained by contracting the tensor 
network is equivalent to the ground state of its corre- 
sponding quantum ID model 18 . The singular point of 
the quasi-entanglement entropy which corresponds to the 
phase transition in quantum ID model can be used to 
determine T c or T c ^tk for the 2D classical lattice. The 
temperature dependence of the quasi-entanglement en- 
tropy Sq at /3 = 1.6 is calculated, as presented in Fig. 
6 (b), where a sharp peak of Sq is seen at the temper- 
ature that is the same as the temperature at which the 
susceptibility exhibits a singularity. It appears that the 
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FIG. 2: (Color online) Tensor network representation in 
LTRG calculations of this present system, (a) the cell lat- 
tice for tensor construction; (b) the schematic representation 
of tensor T; (c) the tensor network of the partition function; 
and (d) the matrix product state (Ref. [12]) obtained by con- 
tracting the tensor network in (c). 
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FIG. 3: (Color online) Temperature dependence of the specific 
heat of the anisotropic 3-state Potts model on the generalized 
Kagome lattice with a = -1. (a) /3 < 1, and (b) j3 > 1. 
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FIG. 4: (Color online) Temperature dependence of the mag- 
netization per site of the anisotropic 3-state Potts model on 
the generalized Kagome lattice with a = -1 and {3 < 1. 
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FIG. 6: (Color online) (a) Correlation length of the 
anisotropic 3-state Potts model on the generalized Kagome 
lattice with a = -1 and /3 > 1, and the inset is for < 1. (b) 
Quasi entanglement entropy and susceptibility of this model 
with = 1.6. 



FIG. 5: (Color online) Temperature dependence of the sus- 
ceptibility of the anisotropic 3-state Potts model on the gen- 
eralized Kagome lattice with a = -1 and f3 > 1. 



quasi-entanglement entropy defined by Eq. (7) can be 
used to detect the KT phase transition. 



IV. REENTRANT PHENOMENA AND PHASE 
DIAGRAM 

The KT phase transition implies the quenching of the 
partial order below T c kt that further supports the dis- 
ordered ground state when (3 > 1 owing to the highly 
degeneracies in the ground state and the strong fluctu- 
ations in the Potts model. Under T Cj kt, the floating 
phase is very sensitive to external perturbations. Fig. 
7 (a) gives the induced m by a tiny field (~ 10~ 3 ) for 
P = 1.6, where one may see that m increases fast with 
increasing the small field, showing the character of the 
floating phase. The inset of Fig. 7 (a) demonstrates that 
m decreases with temperature in a small field, and at a 
certain temperature, has a rapid drop to the value of 5/6, 
indicating again that it is in the floating state. 

From the above discussions, we come to conclude that 
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FIG. 7: (Color online) (a) The magnetization versus small 
magnetic field at temperature T/\Ji\ = 0.32. The inset is 
the magnetization versus temperature in a small field h — 
10 -3 . Here /3 = 1.6. (b) Temperature dependence of m in the 
reentrance region. 



the anisotropic 3-state Potts model on the generalized 
Kagome lattice for (3 > 1 has the topological phase transi- 
tion of the KT type. An interesting question then arises: 
can a reentrant phase exist under a phase without a par- 
tial order through the KT transition? Let us answer this 
question below. 

To search for the reentrant region, we would collect 
the information from the specific heat, susceptibility and 
the quasi-entanglement entropy to determine the phase 
diagrams. We first begin with a = — 1 and take (3 as 
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FIG. 8: (Color online) Temperature dependence of the sus- 
ceptibility of the anisotropic 3-state Potts model on the gen- 
eralized Kagome lattice in the reentrance region. The inset 
shows the corresponding quasi entanglement entropy. 



the variable. In calculations we determine the critical 
temperature of ferromagnetic phase transition by means 
of the peak of the specific heat and the KT phase tran- 
sition temperature according to the peak of the quasi- 
cntanglement entropy. By performing a large-scale cal- 
culations, wc find that the parameter range of j3 where 
the reentrant phenomenon occurs is in between [0.94, 1], 
and in this small region the system is most frustrated. 
Fig. 7(b) gives the temperature dependence of magne- 
tization for a = — 1 and /? = 0.95, and it can be seen 
that with increasing temperature the magnetization de- 
creases sharply to the value of 5/6, after experiencing 
the flat change, then remarkably increases, suggesting 
that there should be two phase transitions in this case, 
namely, one from F state to paramagnetic (P) state, and 
another from P state to floating state. Such transitions 
can also be clearly seen from the susceptibility, as shown 
in Fig. 8, where the temperature dependence of the sus- 
ceptibility is presented for j3 = 0.952 and a = — 1. There 
are two peaks: the first peak at T c j ~ 0.136 corresponds 
to the phase transition from F to P phase, and the sec- 
ond peak at the critical temperature T CtX i m 0.140 in- 
dicates the phase transition from P to X phase. The 
inset of Fig. 8 gives the temperature dependence of the 
quasi-cntanglemcnt entropy, where three peaks are ob- 
served: the two peaks at the critical temperatures 0.136 
and 0.140 are too close to be visibly separated, which just 
correspond to the two peaks of susceptibility, while the 
third peak occurs at T cx f?a 0.242, indicating the phase 
transition from X phase to P phase. Therefore, for /3 
= 0.952, with increasing temperature, the system under- 
goes phase transitions from P to X to P phase, which is 
nothing but the reentrant phenomenon. 

The whole phase diagram for a = —1 is depicted in 
Fig. 9, where three phases including F, P and X phases 
are identified. One may note that there exists a small re- 
gion (3 € [0.94, 1] where the reentrance occurs, as shown 
in the inset of Fig. 9. It can be observed that in this 
region and with increasing temperature, the system goes 
first into the P phase from the F phase, and then enters 
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FIG. 9: (Color online) The phase diagram of the anisotropic 3- 
state Potts model on the generalized Kagome lattice with a = 
-1. The inset indicates the reentrant region. F: ferromagnetic 
phase; P: paramagnetic phase; X: Floating phase. 
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FIG. 10: (Color online) Phase diagram of the anisotropic 3- 
state Potts model on the generalized Kagome lattice with < 
a < 1. F: Ferromagnetic; P: Paramagnetic; X and Y: floating 
phases. The inset is a-dependence of the ferromagnetic phase 
transition critical temperature T c j. 



into the X phase, and then goes into again the P phase. 
In addition, some further computations reveal that the 
region of (3 for which the reentrant phenomenon can ap- 
pear is proportional to |a| (a < 0) and shrinks to zero at 
a = that is just the typical Kagome lattice. It is also 
worth mentioning that the phase diagram in Fig. 9 is 
similar to that of the generalized Kagome Ising model 3 . 
Since in the 3-state Potts model the thermal fluctuations 
are stronger than in the Ising model, the phase bound- 
aries in Fig. 9 stand lower than in the diagram of Ref. 
[3]. In the latter diagram the region of the partial or- 
dered state A is replaced by the floating phase labeled 
by*. 

When < a < 1 , there is also a reentrant property, as 
shown in Fig. 10. In the phase diagram with a = 0.25, 
there are four phases, including F, P, X and Y phases, 
where Y is also a floating phase appearing in a very small 
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FIG. 11: (Color online) Phase diagram of the anisotropic 3- 
state Potts model on the generalized Kagome lattice with a 
> 1. P: Paramagnetic; X and Y: floating phases. 



region when ft > 0.8. There exists a reentrant region 
between F and X phases when ft < -24.94. This region, 
determined by the same procedure as above, is the same 
as the generalized Kagome Ising model 3 ' 5 , but the critical 
temperature T c j is lowered than that in the Ising case 
due to strong thermal fluctuations. It is found that the 
area of X phase increases with the value of a in the cost 
of the decreased scope of the F phase. It is clear that in 
the inset of Fig. 10 the ferromagnetic phase transition 
temperature T c j goes to zero when a increases near to 
1. The finite temperature phase transition from X to P 
phase is also of the KT type, and the critical temperature 
is obtained through the calculation of the entanglement 
entropy or susceptibility. 

When a > 1 the F phase disappears, and there will be 
only two floating phases X and Y in the phase diagram, 
as shown in Fig. 11. In this case, no reentrant phe- 
nomenon occurs. In addition, the partial ordered state 
with configuration B in the Ising model (Ref. 3) when a, 
ft > is vanished due to highly thermal fluctuations in 
the 3-state Potts model as in Fig. 10. Besides, we found 
through a number of calculations that when q > 3, the 
g-state Potts model on the generalized Kagome lattice 



exhibits no any reentrant property. 



V. SUMMARY 

To summarize, we have studied the anisotropic 3-state 
Potts model on the generalized Kagome lattice by utiliz- 
ing the LTRG method. The phase diagrams in the plane 
of temperature versus ft for the cases ofa<0, 0<a<l 
and a > 1 are obtained. Different phases including fer- 
romagnetic, paramagnetic and floating phases are identi- 
fied. The phase boundaries are determined by observing 
the singularities in the specific heat, magnetization, sus- 
ceptibility and quasi-entanglement entropy. For the two 
cases of a < and < a < 1, the reentrant phenomena 
are observed in small regions, respectively. By studying 
the behaviors of the correlation length, we found that 
the phase transition for ft > 1 is of the KT type, in sharp 
contrast to the case of ft < 1 where the correlation length 
has a quite different behavior from that of ft > 1. 

By comparing with the generalized Kagome Ising 
model, where the partial ordered phases exist, the present 
3-state Potts model possesses floating phases X and 
Y owing to strong thermal fluctuations. Such floating 
phases do not have any local order parameters. Through 
the present study, we may remark that the appearance of 
the reentrant phenomena in 2D classical lattices is mainly 
caused by frustrations as well as the geometrical struc- 
ture of the lattices, which has also close relation with the 
freedom q of the local spin. Importantly, the existence of 
the partial ordered state in the ground state may not be 
the necessary condition for the reentrant phenomena. 
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